The intriguing life of massive galaxies 
Proceedings IAU Symposium No. 295, 2013 
D. Thomas, A. Pasguali & I. Ferreras, eds. 



©2013 International Astronomical Union 
DOI: 00.0000/X000000000000000X 



Modelling the formation of today's massive 

ellipticals 

Thorsten Naab 1 

1 Max-Planck-Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 85741 Garching, Germany 

email: naab@mpa-gar ching . mpg . de 

Abstract. The discovery of a population of massive, compact and quiescent early-type galaxies 
has changed the view on plausible formation scenarios for the present day population of ellip- 
tical galaxies. Traditionally assumed formation histories dominated by 'single events' like early 
collapse or major mergers appear to be incomplete and have to be embedded in the context of 
hierarchical cosmological models with continuous gas accretion and the merging of small stellar 
systems (minor mergers). Once these processes are consistently taken into account the hierar- 
chical models favor a two-phase assembly process and are in much better shape to capture the 
observed trends. We review some aspects of recent progress in the field. 
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1. Introduction 

During the formation and assembly of massive galaxies merging is a natural process 
in modern hierarchical cosmological models. It is expected to play a significant role for 
the structural and morphological evolution (e.g. Kauffmann et al.lll996t Kauffmannlll996 



De Lucia et al.l2006tlKhochfar k Sil k 2006; De Lucia k Blaizotl2007l:lGuo k WhiteT2 008: 



Kormendv et al.l l2009t Hopkins et al.l l2010af) . In the light of these theoretical expec- 



tations and d irect observations of 'dry' mergers of gas poor ellipt ic al galaxies up to 
high redshift (Ivan Dokkurnll2005t iTran . et al.| [2005l: iBell et al.ll2006al lbl; lLotzll2008t Ijogee 
20091: Newman et al.ll2012t iMan et al.l l2012h simulations of idealized collisionless merg- 



ers have again received attention and new studies were triggered. Merger simulations 
of already existing spheroidal galaxies have focused in detail on the evolution of abun- 
dance gradients, shapes and kinematics, s caling relations, sizes and dark matter fractions 
( White! Il978l 119791: iMakino k HutllT997t iBovlan-Kolchin et al.ll2005t iNaab et al.ll2006b : 
Bovlan-Kolchin et al.l 120061 12008: Pi Matteo et al.l 120091 : iNipoti et al.l " l2009bl l2012aF 7lf 



the progenitors were two disk galaxies (which then can include a gaseous component) 
the aim was to investigate the morphological transformation, i.e. the formation of new 
dynamically hot spheroidal elliptical galaxies from two dynamically cold progen i tor spi- 
ral galaxies (lGerhardlll98lt iFarouki k Shapirolll982t iNegroponte k Whitej|l983l; iBarnesI 
1 19881: iBarnes k Hernauist 1992t Hernauist 1992). Apart from studies of the effect of 
the merger mass-ratio (lBarneslll998l:lBekkilll998l:lBendo k Barnesll200fl iNaab k Burkert 
l2003HBournaud et al.ll2004 120051: iGonzalez-Garcia k Balcellsll2005l ) the tidal torquing of 
gas, its inflow to the central regions, the impa ct on the stellar orbits (IBarnes k Hernauist 
1996 | :lNaab et a l. 2006a; Ho ffman et al.ll2010f). su bsequent starbursts ( Mihos k Hernauist 



199411 996: Barnes 2004; Di Mat teo et al.ll2008l) and the potential growth of black holes 



(lHernauistjl989tlSpringel et al.l2005l:lDi Matteo et al.ll2005l I Johansson et al.l2009aHYounger et al 



2009) was investigated in numerous studies together with influential studies on the 



origin of early-type galaxy scaling relation s (iRobertson et al 



Cox et al.ll2006l lHopkins et al.ll2008l l2009dH lai nDebuhr et al 



2006; Dekel k Coxi 2006 



20101: iMoster et al.ll201ll ) 
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However, despite the detailed insights on the stellar and gas dynamical processes in sim- 
ulated galaxy mergers, the 'binary merger' approach is limited in scope and seems not 
to be able to naturally explain all properties of present day massive elliptical galaxies 
( Naab fc Ostrikerll2009l ). 

The most massive elliptical galaxies (or their progenitors) are considered to start 
forming their stars at high redshift (z ~ 6, or hig her) in a dissipative environment 



Feldmann et al.l 120 lit lOser et al 



20121) . 



rapidly become very massive (~ 10 L1 M^,) by z = 2 ( Keres et al. 20051 Khochfar fe Silk 
20061 iDe Lucia et al.ll20o"i iKriek et alJbooa [Naab et al.ll2007l. l2009t Ijoung et alJl200a 
Dekel et al.l2009HKeres et al.|2009HOser et al.l2O10HFeldmann et al.ll2010l : lDoim'nguez Sanchez! 



A significant fraction of this high red- 
shift population is observed to be already quiescent at z ~ 2, on average 4-5 times 
more compact (part o f this apparent evolution might driven by selection effects, see e.g. 



Poggianti et al.ll2012f). and typically a factor of two less massive than their low redshift 
descendants (Daddi et al.l 1200 



Truiillo 2006; Longhetti et al 
20081 Ivan der Wei et al ~ 



io5 Ivan der Wei et alJEioll Idi Serego Alighieri et al.l 12005 
. 2007tlToft et alJ2007tlBuitrago et a .l2008Hvan Dokkum et 



Buitrago 

2008t iCimatti et al.l [2008 ; jFranx et al.l l2008t iDamianov et al. 
200S ; ICenarro fc Truiilloll2009l;lBezanson et al.lbooalvan Dokkum et al l2O10Hvan de Sande eTal 



et al 



201lHWhitaker et al.ll2012h . It is reasonable to assume that the high-redshift population 
forms the cores of at least some, if not all, present day massive ellipticals. This rapid 
structural evolution is supposed to happen in an inside-out fashion, mainly by adding 
stellar mass to the outer parts of the ga laxies over time, however, without the forma- 
tion of a significant fraction of new stars (|Hopkins et a l. 2009a; Ivan Dokkum et"alll2010t 
Szomoru et al. 2012; Sa racco et al. 2012). In this respect the growth of massive quiescent 
high-redshift galaxies is markedly different to the star formation driven inside-out growth 
of disk galaxies 

The implications of these observational findings for the formation and evolution of 
massive elliptical galaxies are many-fold. They are unlikely to have formed by an initial 
'monolithic collapse' followe d by passive evolution as their present day counterparts would 
be too small and too red (Ivan Dokkum et al.l 12008 ; Kriek et al. 20081 : Bezanson et al 



20091 : Ferre-Mateu et al~ll2012 ). In addition the evolution of these system cannot be ex- 
plained by just a single 'binary merger of disk galaxies'. The compact high-redshif t sys- 
tems might have formed in such a process ( Wuvts et al. 20101 Bournaud et al. 2011 ). if it 
were gas-rich, but the subsequent structural evolution requires additional processes which 
are not driven by the formation of new stars. Observational results that a lmost none of 



these massive compact galaxies were able to survive to the present day (|Trujillo et al. 
20091: iTavlor et al.ll2010h indicate that a general and common physical mechanism must 



be at work. Spectacular events alone, like major early-type galaxy mergers, might be too 
rare. 



2. Minor mergers vs. major mergers 

Minor merges, however, are expected to happen frequently in the lifetime of a massive 
galaxy and have received particular attention as they provide a natural way to increase 
the size of a galaxy. With only a few assumptions the virial theorem provides a sim- 
ple estimate of how a one-component system evolyes du ring major and minor mergers 
(|Cole et alJboOOl : iNaab et al.l 120091 iBezanson et"~aT1l2009h . Following iNaab et al.l (|2009l ) 
we assume that a compact initial stellar system has formed (e.g. involving gas dissipa- 
tion) with a total energy a mass M^, a gravitational ra dius r a .j, and the mean sq uare 
speed of the stars is (vf ). According to the virial theorem (jBinnev fc Tremainell2008l) the 
total energy of the system is 
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Figure 1. Left: Simulated size evolution as a function of bound stellar mass for mergers with 
mass-ratios 1:1 (blue), 5:1 (red), and 10:1 (green) . The observationally expected relation is 
indicated by the black line ijvan Dokkum et afl [2010h The presence of dark matter significantly 
boosts the size evolution of 5:1 and 10:1 mergers. Right: This also lea ds to a significantly stronger 
evolution of the Sersic index (figures taken from lHilz et al.l (|2012aT 



Ei = Ki + Wi = -Ki = 

1 . o, 1 GMf 

This system then merges (on zero energy orbits) with other systems of a total energy 
E a , total mass M a , gravitational radii r a and mean square speeds averaging (v%). The 
fractional mass increase from all the merged galaxies is rj = M a /Mi and the total kinetic 
energy of the material is K a — (l/2)M a (v^) , f urther defining e = (vjj)/(vj ). Under the 
assumption of energy conservation (results from ( Khochfar fc Burkertll2006l ) indicate that 



most halos merge on parabolic orbits) the ratio of initi al to final mean square speeds, 
gravitational radii and densities can be then written as (|Naab et al.ll2009h 



(vj) _ (l + r,e) r gJ _ (1 + V f p f _ (1 + V ef 



(vf) l + r/ ' r g>i (1 + r?e) ' (1 + n) 5 ' 

For mergers of two identical systems, rj = 1, the mean square speed would remain 
unchanged, the size increases by a factor of two and the densities drop by a factor of 
four. In the limit that the mass is accreted in the form of very weakly bound stellar 
systems with (v^) << (vf) or e << 1, the mean square speed is reduced by a factor 
two, the size increases by a factor four and the density drops by a factor of 32. These 
estimates are, however, idealized assuming one-component systems, no violent relaxation 
a nd zero-ene r gy orb its with fixed angular momentum. 



Hilz et al.l (|2012bl ) have recently re-investigated in detail the collisionless dynamics of 
major and minor mergers of systems including concentrated stellar spheroidal compo- 
nents embedded in extended dark matter halos. They present more accurate versions of 
the above equations including the effect of escapers and the interaction of the stellar bary- 
onic with dark matter and describe in detail how the presence of a massive dark matter 
halos alter the evolution of the merging systems. One result of this study was that both 
minor and major mergers lead to size growth and an increase of the dark matter fraction. 
The physical processes are, however, different. Violent relaxation in major mergers mixes 
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Figure 2. Left: Size evolution driven by rapid mass-loss from the idealized simulations of isolated 
galaxies. If 40 per cent of the mass is lost (e = 0.6) the sizes can rapidly increase by 60 per cent. 
From the black to th e green line the ejection varies from i mmediate ejection to an ejection time of 
80 Myrs (taken fron yRagone-Fi g ueroa &i Gra nato 2011). Right: Evolution of the stellar surface 
density profiles of a cosmological zoom-simulation of a brightest cluster galaxy in a model with 
strong AGN feedback. Due to the gas explusion from the AGN the system is sig nificantly more 
exten ded than in the no- AGN case and even develops a central core (taken from|Martizzi et al.l 
12012ft . 



dark matter to the central regions. Escaping, unbound, particles limit the expected size 
growth to values below the ones expected from the idealized equations above. In minor 
mergers (mass-ratios of 1:5 and 1:10), the stellar satellites are stripped at large radii 
where the host galaxies dominated by dark matter and the stellar effective radii and the 
dark mat ter fractions grow m ore rapidly than expected from the simple virial equations 
(see a lso Laporte et al.ll2012 ). Due to the addition of stellar satellite material at large 



radii ( Villumsenl 1983 ). the stellar mass distribution changes significan tly resulting in a 
significant increase of the Sersic index (see Fig. [T1 and iHilz et al.ll2012al). The general re - 
sults on size evolution are in agreement with similar studies by e.g. lOogi k. Habd ( 20121 ). 
However, there is an ongoing debate of whether the size growth by minor mergers is suf- 
ficient to explain th e observed cosmological size evolution of elliptical galaxies. Whereas 



Oogi fc Habel d2012h argue th a t the s iz e growth by minor m erger s alone might be suffi - 



cient, studies bv lNipoti et al.l (|2012b[ ). ICimatti et al.l (|2012l ). and iNewman et al.l (|2012f ) 
have combined idealized numerical simulations embedded in a cosmological context and 
new observational constraints. They come to the conclusion that minor mergers might 
be able to explain the observed size growth from redshift z ~ 1 to the present. However, 
at higher redshift minor and major mergers might not be frequent enough to explain the 
rapid size evolution observed at z > 1 and therefore an additional physical mechanism 
might be required. 

A potential candidate for such an additional process is AGN driven outflow of gas 
from a massive high-re dshift gas-rich and compact galaxy ( Fan et al1l2008t Hopkins et al 



2010btlFan et alfc oiO). In general, ste llar systems suffering from central mass-loss ei 



-Mfinai /-Winitiai wm expand (Hills 198011 and for rapid and slow mass-loss simple relations 
for the ratio of the final to the initial radius can be derived: 
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Figure 3. Upper left: Examples for the assembly history (stellar origin) of three massive galaxies 
in high-resolution cosmological zoom simulations. At high redshift the formation is dominated by 
in-situ star formation (red colors ). The low redshift ass embly is dominated by merging of stellar 
systems (blue colors, taken from iFeldmann et al.l l2010). Upper right: Ratio of the accreted over 
final stellar mass versus final stellar mass for galaxies in a cosmological simulation box (void: 
blue, cluster: red) including strong supernova feedback (upper panel). The fraction of accreted 
stars is about a factor 2 -3 lower than in the high- resolution zoom simulations of lOser et al.l 
(|2010D without strong supernova feedback (lower panel); the trend with mass is similar but 
less strong (taken from lLackner et al.ll2012T h Lower left: Independent estimate of the ratio of 
accreted to in-situ formed stars as a function of halo mass from abun dance matching studie s 
ijMoster et al.|[2012f ). Lower right: Similar estimates from a study by iBehroozi et all (|2012f ). 
Both studies find a strong trend that the assembly of galaxies in more massive halos is more 
dominated by the accretion of stars rather than in-situ star formation. 



Rfinal,rapid _ ^lo 



^initial 2q oss 1 
Rfinal,slow 1 



toss 



Rinitial £lo 

It is worth noting that rapid mass-loss of more than half the total mass can un- 



6 



Thorsten Naab 





Figure 4. Upper panels: The present day mass-size relation (left) for a sample of high-resolution 
zoom simulations (blue points, full symbols are quiescent galaxies) compared to observations. 
The evolution of the relation is driven by accretion of stars and is indicated by the location of 
the most massive progenitor galaxies at different redshifts. For all galaxies more massive than 
the mass limit indicated on the left plot (log(M„) = 10.8) the average size evolution agrees well 
with observations. Lower panels: Similar plot for the evolution of the mass-dispersion relation 
(left). In a fixed mass range galaxies have higher dispersions at higher redshifts (rig ht). Again, 
the s imulated evolution is very similar to the observed one (figures are taken from lOser et aU 

Eo]J). 



bind the whole system. Thi s process is well known and has been studied for star clus- 
ters (lHillslll980l). galaxies (lHilldll980HHopkins et^l2010bl;[Ragone-Figueroa fc Granato 



2011 



Pontzen fe Governato 2012I ). as well as cores of galaxy clusters (see Fig. f2] and 



Martizzi et al.ll2012ft . 



3. The cosmological two-phase assembly 

The assembly histories of massive galaxies in currently favored hierarchical cosmolog- 
ical models are significantly more complex than a single binary merger. They grow - in 
particular at high redshift - by smooth accretion of gas, major mergers but also numerous 
minor mergers covering a large range of mass-ratios which can dominate the amount of 
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assembled stars. The picture that is emerging from from semi-analytical models and high- 
resolution cosmological simulations of massive galaxies bears a two-phase characteristic 
(iDe Lucia fc Blaizodl2007t iGuo fc White! 120081: iGenel et alj|2008t iFeld mann et alj |20ld: 
Oser et al.ll2010HFeldmann et alj|201lb iHirschmann et al. 20121) . 

At high redshifts the formation is dominated by dissipative processes (i.e. significant 
radiative energy losses) and in-situ star formation leading compact progenitors with high 
phase space densities. In a second phase massive galaxies are growing by the addition 
of stars at large radii that have formed early outside the main galaxies in other galaxies 
that were accreted later-on. This assembly phase is dominated by collisionless dynamics 



and radiative energy los s es are of minor imp ortance (see e.g. iJohansson et al 
Lackner fc Ostrikeill2oTot lLaporte et alj|2012h . 



2009b 



Independent studies using cosmological simulations based on different numerical meth- 
ods come to similar conclusions that -on average - the m ass assembly of massive galaxies 



is dominated by minor mer gers with mass-ratios ~ 1 : 5 (jOser et alj|2012t lLackner et al 



20121: iGabor fc Pavel 120121 ) The relative importance of accreted versus in-situ formed 
stars inc reases with galaxy mass, a result that was already predicted by semi- analytical 
models (|De Lucia et al.l 120061: IDe Lucia fc Blaizotl l2007t IGuo fc White] |200j) and has 
been confirmed by indepen dent estimates from abundance matching techniques ([Moster et al 
2012HBehroozi et al.ll2012l) . The absolute fractions are model dependent and can vary e.g. 
by ~ 50% for different feedback models (see Fig. [3j) . Studies based on cosmological zoom 
simulations make a plausible point that the present day scaling relations might be set 
by the stellar accretion history of massiv e galaxies, i.e. the above mentioned fraction of 
in-situ to accreted stars (Os er et al 1 l2012h . In addition, based on still small samples, high- 
resolution cosmological simulat ions the evolution of the scaling relation s appears to be in 



accordance with observations (Feldmann et al. 2011; Oser et al. 2012; Johansson et al 



20121 ). However, in general the cosmological simulations of massive galaxies still fail to 
reproduce all observational constraints at the same time and are still limited with respect 
to either resolution and statistics as well as the algorithmic implementation of relevant 
feedback processes. In particular feedback from super-massi ve black holes might help to 
finally meet observational c onstraints for massive ellipticals (jMcCarthv et al.ll20ldl20lil 
Puchwein fc Springelllioij ). 



TN acknowledges support by and valuable disucssions with Peter Johansson, Ludwig 
Oser and Jeremiah P. Ostriker. 
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